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Saccharomyces cerevisiae possesses numerous advantageous biological features, such as being robust, easily handled, 
mostly non-pathogenic and having high catabolic rates, etc., which can be considered as merits for being used as a 
promising biocatalyst in microbial fuel cells (MFCs) for electricity generation. Previous studies have developed efficient 
MFC configurations to convert metabolic electron shuttles, such as cytoplasmic NADH, into usable electric current. 
However, no studies have elucidated the maximum potential of 5. cerevisiae for current output and the underlying 
metabolic pathways, resulting from the interaction of thousands of reactions inside the cell during MFC operation. To 
address these two key issues, this study used in silico metabolic engineering techniques, flux balance analysis (FBA), 
and flux variability analysis with target flux minimization (FATMIN), to model the metabolic perturbation of S. cerevisiae 
under the MFC-energy extraction. The FBA results showed that, in the cytoplasmic NADH-dependent mediated electron 
transfer (MET) mode, S. cerevisiae had a potential to produce currents at up to 5.781 A/gDW for the anaerobic and 6.193 
A/gDW for the aerobic environments. The FATMIN results showed that the aerobic and anaerobic metabolisms are 
resilient, relying on six and five contributing reactions respectively for high current production. Two reactions, catalyzed 
by glutamate dehydrogenase (NAD) (EC 1.4.1.3) and methylene tetrahydrofolate dehydrogenase (NAD) (EC 1.5.1.5), were 
shared in both current-production modes and contributed to over 80% of the identified maximum current outputs. It is 
also shown that the NADH regeneration was much less energy costly than biomass production rate. Taken together, our 
finding suggests that S. cerevisiae should receive more research effort for MFC electricity production. 



Introduction 

Saccharomyces cerevisiae is associated with anthropic environ- 
ments and is well known for food or beverage (alcoholic) fer- 
mentation. 1 Besides these applications, this yeast species has 
been employed as an eukaryotic model organism in molecu- 
lar and cell biology; for example, the characteristics of many 
proteins can be discovered by studying their homologs in 
S. cerevisiae. 2 ' 6 Recently, it has been shown that S. cerevisiae 
can be processed to produce potential advanced biofuels such 
as long chain alcohols and isoprenoid- and fatty acid-based 
biofuels, which have physical properties that more closely 
resemble petroleum-derived fuels. 7 Nevertheless, biofuels need 
to be further combusted to produce usable energy in the form 
of electricity. To circumvent the disadvantages of biofuel com- 
bustion and directly convert the metabolic reducing potential 
inside the cell into electricity, another device, named a micro- 
bial fuel cell (MFC), has also been proposed. In a typical MFC 
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configuration, microorganisms are grown in the anodic com- 
partment and obtain electrons from breaking down substrate 
during growth. 8 The gained electrons can be harnessed and 
transported out of cells to an electrode via two main opera- 
tional modes, mediated electron transfer (MET), and direct 
electron transfer (DET) modes. 8 MET mode involves soluble 
redox mediators (such as bromocresol green [BG] and neutral 
red [NR] 9 ) that act as electron relays repeatedly cycling between 
the metabolism of the cell and the electrode, 10 whereas DET 
mode is the case where electrons are directly expelled by some 
microorganisms via membrane-associated proteins (e.g., c-type 
cytochrome). 10 In either of these modes, the physiological effect 
of the immersion of a microbe in the anodic compartment is 
that an electron drain is introduced by the external environ- 
ment. Because of the key role of redox carriers in most aspects of 
metabolism, operation of the metabolic network is fundamen- 
tally affected by this. The focus of this article is to model the 
resulting rebalancing of metabolic fluxes by using flux balance 
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analysis (FBA). Further exposition of the working principle of 
an MFC can be found in several review articles. 8,11 ' 12 

Both MET and DET modes can be used in MFCs based on 
S. cerevisiae." In the DET mode, the yeast based MFC could so 
far only produce an extremely low current output, at a scale of 
microamperes. 13 In addition, it is unclear what mechanisms on 
the cell wall are responsible for DET current production. 14 In the 
MET mode, mediators are required to facilitate the transfer of 
electrons to the anode, because this yeast is thought incapable 
of producing such mediators indigenously. 15 Because glycolysis 
occurs in the cytosol of the yeast cell rather than in the mito- 
chondria, the electron shuttle, NADH, is easily accessible to a 
mediator molecule present in the cell membrane of the yeast and 
has been proposed as the main electron supplier in a yeast-based 
MFC. 15,16 Based on these characteristics, MFCs using yeast can 
be directly applied in fermenters for in situ power generation. 15 
Other electron carriers (such as FADH,) may also contribute to 
current production, but these reduced molecules are regenerated 
by much fewer reactions than NADH, which is widely used in 
many reactions of the energy metabolism (e.g., glycolysis and 
Krebs cycles). 

Since the biological entity, S. cerevisiae, in the MFCs are con- 
sidered as the fundamental part that determines the current and 
power output of MFC and draws a distinction between MFC and 
other chemical fuel cells, it would be desirable to know the inher- 
ent restrictions and mechanisms of the yeast for current genera- 
tion and to understand how the energy extraction perturbs the 
metabolic mechanisms of microorganisms in MFCs. Although 
some studies have discussed the metabolisms for supplying elec- 
trons in MFCs based on empirical knowledge of the textbook 
biochemical pathways or limited results arising from reductionist 
investigative methods, there is still an urgent demand for knowl- 
edge about the metabolic pathways resulting from the interaction 
of thousands of reactions inside the yeast cell during MFC opera- 
tion. This task can be solved by in silico modeling of genome-scale 
metabolic network reconstructions (GENRE), using the recently 
developed constraint-based modeling approaches. 17 ' 18 However, 
these novel methods had not been used in MFC research area 
until recently. 11 ' 19,20 

In this study, to identify the maximum potential of S. cere- 
visiae to sustain the energy extraction process in yeast-based 
MFCs, FBA, and multi-objective optimization 19,20 was employed 
to investigate the metabolic capability of the yeast to supply an 
excess flux of cytosolic NADH to reduce the oxidized media- 
tors (or anode), under two broad metabolic types, anaerobic and 
aerobic growth. A genome-scale metabolic model of S. cerevisiae 
(yeast 5) was chosen as the backbone for the present flux balance 
modeling, since it has been tested extensively against experimen- 
tal results for prediction of several growth conditions, including 
byproduct synthesis, in a previous study. 21 Only minimal changes 
were made to the model for our study, to represent the diversion 
of redox carriers in a MFC. In addition, flux variability analysis 
with target flux minimization (FATMIN) 20 was implemented to 
elucidate pertinent reactions underlying the maximum current 
output. Finally, robustness analysis was conducted to establish 



the feedstock costs for the cytosolic NADH-linked electricity 
generation in S. cerevisiae. 

Results and Discussion 

Impact of the redox perturbation on the biomass production. 
Figure 1 shows how production of mediated electron transfer 
(MET) aimed at cytosolic NADH competes with biomass pro- 
duction for metabolic resources. MET influenced the production 
of the NADH and biomass in the same way for both anaero- 
bic and aerobic conditions, that is, the increase in MET drove 
the NADH regeneration rates in the two growth modes toward 
their maximum attainable values, and the corresponding NADH 
flux consumed for cell maintenance and biomass formation rates 
toward zero. This indicates most of the energy originally for 
growth is converted to electricity. 

At the base state optimized for growth, the total turnover rate 
of NADH were 19.74 mmol/gDW/h and 67.34 mmol/gDW/h 
for the anaerobic and aerobic condition respectively. Compared 
with the base states, the anaerobic and aerobic metabolisms had 
a potential to increase their NADH regeneration rates by about 
5.469 (447%) and 1.716 (71.6%) folds respectively under the 
highly NADH-perturbed metabolic states, achieving the high- 
est NADH_mfc flux values of about 108 mmol gDW" 1 h" 1 for 
the anaerobic growth and 115.6 mmol gDW" 1 h" 1 for the aero- 
bic growth, which were achievable when biomass production 
rate was suppressed to 0.012 and 0.01284 h" 1 in the anaerobic 
and aerobic modes respectively. This indicates that aerobic res- 
piration can increase the maximum biomass production rate, 
but does not substantially influence the metabolic capability 
for NADH regeneration. At the same substrate consumption 
rate, the fermentation can regenerate NADH to a level that is 
sufficient to sustain similar electric current output as aerobic 
growth. 

The MET production rates shown in the top parts of Figure 1 
for the two growth cases, only show a linear dependence on bio- 
mass growth for small growth rates. Starting from the unper- 
turbed state of zero MET production at the right hand side of 
the plots and based on the slope changes in the lines connecting 
two neighboring data points, the relationships between varying 
biomass production rate and the MET (NADH_mfc) rate can 
be divided into three stages (see Fig. 1). These are: an increase in 
MET without changing biomass growth (Stage I), a downward 
curving nonlinear rise of MET vs. growth rates (Stage II), and a 
linear dependence of the rise in MET on drop in biomass forma- 
tion (Stage III). Stage III only appears in the left corner of the 
figures, where the line connecting the points clustered indicates 
there is a linear relationship between the biomass production and 
the current output. 

The slope of the NADH production vs. biomass growth curve 
as in Figure 1 gives a direct measure of the efficiency of the 
underlying pathways. A large slope means that NADH produc- 
tion increases by a large amount for a given cost in terms of bio- 
mass production and is therefore highly efficient and vice versa. 
In the discussions below, we base interpretations of reallocations 
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Figure I.The relationships of the biomass production rate and cytosolic NADH flux diverted toward anode (NADH_mfc). The NADH supplying rate 
in the MET mode and the reducing equvalent consumption rate for cellular use, as functions of biomass production rate. The red line represents the 
maximal NADH_mfc flux for a feasible biomass production rate, while any point within the pink area represents all allowable NADH_mfc and biomass 
production rates. The blue area represents the total cytosolic NADH consuming flux for normal cellular function. The distance between the two lines 
across the pink and blue areas represents the total available reducing equivalent flux in the cell at a metabolic state related to a specific biomass pro- 
duction rate; inset, enlargement of boxed area. Dashed line with arrowheads indicates which two data points are paired up. The points are simulated 
by bi-objective optimization involving varying coefficients (X) assigned for the growth and MET maximizations. A more detailed discussion of the 
significance of the X is given in the Methods section. Briefly, it is a parameter that continuously adjusts the metabolic state of the cell through a range 
stretching from pure growth without the extraction of current (at one extreme) to at state in which all metabolic resources are taken up by supplying 
electric current and no growth, at the other extreme. 



of metabolic resources to pathways with different efficiencies on 
this evidence directly derived from the modeling results. 

Looking at the net NADH production (MET) curve, the 
Stage I relationship indicates that solutions that are based on 



reallocation of energy resources inside cells can elevate the slight 
increase in MET rate, which obviate the degradation of the 
intrinsic biological objective, to grow. However, rerouting the 
energy flow through different metabolic pathways could quell 
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some cellular functions, which, for example, are unnecessary in 
optimal growth conditions, but may be important to survival 
in adverse conditions. This also implies that maximization of 
growth rate is not the sole goal for eukaryote. 

The Stage II indicates that the metabolic pathways have dif- 
ferent efficiencies and capability in NADH regeneration. Some 
pathways possess higher efficiencies in converting the metabolic 
resource into NADH but lower maximum capability in turnover 
rate and some others vice versa. When the current output was 
relatively low, the routes with higher efficiencies were selected 
to minimize the metabolic resource cost in order to maximize 
growth rate. When the current output rose to a higher level, those 
pathways with lower efficiencies but higher upper bounds of the 
NADH flux had to be used instead. 

The Stage III behavior begins at very low growth rate values 
(less than about 10% of the highest growth rates) in each of the 
two modes. The Stage III results from the fact that, only a certain 
set of pathways are able to sustain the very high current output 
and the flux through these pathways are linearly proportional to 
current output and biomass production rate. 

Figure 1 also shows that a rising MET production rate reduce 
the cytosolic NADH consumption rate until it finally reaches 
zero, which indicates that no cytosolic NADH fluxes were con- 
sumed for cell maintenance. The curve plotted for the NADH 
consumption rate vs. the biomass production can be classified 
into three types: Type I is a reduction in the NADH consump- 
tion rate without changing the biomass production rate; Type II 
represents a decrease in the NADH consumed for maintenance, 
associated with a drop in biomass production rate; Type III is 
the phenomenon that (nearly) no cytosolic NADH flux was con- 
sumed for cell maintenance accompanied by a decline in the bio- 
mass production. 

A correspondence can be found between the types of con- 
sumption behaviors and the stages displayed by the MET pro- 
duction curves. Stage I coincides with Type I, which is reflected 
by the observation that the increase in the MET is accompanied 
by a drop in the NADH consumption for cellular maintenance. 
This indicates that decreasing NADH flux for internal use would 
be the primary mechanism implemented by the yeast to relieve 
the redox perturbation. The similar mechanism has also been 
seen for the NADH-targeted MET mode of C, reinhardtii, 21 This 
suggests that excess cytosolic NADH may be a common occur- 
rence in eukaryotic microorganisms. 

The Stage II occurs with both Type II and III. This indicates 
that in Stage II, there is a continuous adjustment in the metabolic 
strategies underlying the MET current production. The reduc- 
tion in NADH consumption for maintenance may contribute to 
the major portion of the elevated NADH flux diverted to the 
mediator in MFCs in the beginning (i.e., the metabolic states 
undergoing a relatively lower degree of the current-producing 
perturbation). As the NADH flux shunt further augments, the 
intracellular pathways that consume cytosolic NADH are disfa- 
vored by the metabolism. The cell opts for metabolic pathways 
that can excessively produce NADH flux for oxidized media- 
tors, while trying to maximize biomass production rate. When 
entering Stage HI of the MET, the cell redirects nearly all of the 



metabolic resource to NADH regeneration pathways that can 
supply a high flux of NADH. 

Both Stage I and Type I behaviors occurred in the cases of 
C. reinhardtii 12 and S. cerevisiae. The coincidence indicates 
that the multiplicity of eukaryotic functions make these two 
microorganisms tolerable to hostile situations, since it is origi- 
nally reflected that there are more routes of energy flow under- 
lying more physiological conditions. In addition, different from 
C. reinhardtii, where the shape of cytosolic NADH consumption 
line changed in a complex manner, the two cases of S. cerevisiae 
is much simpler. This discrepancy may be a result of the photo- 
synthesis which contains a number of energy pathways providing 
metabolic fuel to the cell growth. 

To further elucidate which one of the two growth conditions 
can better cope with the NADH-related perturbation, we per- 
formed fractional benefit analysis, which calculated a measure 
called the fractional benefit B as discussed in the methods sec- 
tion. Briefly, B is the sum of the fractions of maximal growth and 
NADH production respectively, achieved in a particular meta- 
bolic state. Figure 2 shows the result of this measure applied to 
the reported simulations. Starting from wild type growth rates 
(right hand side of the figures), an increase in B values can be 
seen for the Stage I and Type I behavior. The apex points of the 
B value curves correspond to the highest achievable combined 
benefit. The maximum B value for the anaerobic growth is a 
little bit higher than that for the aerobic growth. This suggests 
that, even though the aerobic growth has a potential to funnel 
more excess NADH for current production than the anaero- 
bic condition, but the loss of NADH can cause more adverse 
effects on the aerobic than the anaerobic growth. Both growth 
conditions can yield a maximum B value of above 75%, which 
is much higher than those of all other three microorganisms 
(G. sulfurreducens, C. reinhardtii, and Synecbocystk sp. PCC 
6803) that we studied previously. 19 ' 20,22 This indicates that addi- 
tional NADH perturbations would result in less hostile effects 
on S. cerevisiae than the other three microorganisms and thus 
S. cerevisiae has the greatest suitability for current production 
among the four. 

The metabolic strategies supporting the high NADH regen- 
eration rate in the anaerobic and aerobic modes are elucidated in 
the next section. 

Metabolic strategies for increasing flux of reducing equiva- 
lents. A list of identified reactions (enzymes) that were respon- 
sible for the promoted diversion of the cytosolic NADH flux 
toward current production is summarized in Table 1. 

Anaerobic growth. In the anaerobic fermentation, S. cerevi- 
siae relied on six reactions to regenerate NADH from NAD* at 
the high rate, subject to the stoichiometry and biological con- 
straints applied. Two of six NADH production routes, namely 
glutamate dehydrogenase (NAD) (EC: 1.4.1.3) and methylene- 
tetrahydrofolate dehydrogenase (NAD) (EC: 1.5.1.5), were the 
main NADH suppliers. The two enzyme catalyzed reactions 
were capable of producing up to 77.48% (85.37 mmol/gDW/h) 
and 81.23% (85.37 mmol/gDW/h) respectively of the maximum 
net NADH supplying rate (105.1 mmol/gDW/h). Nevertheless, 
only the combined fluxes of these two reactions that made up any 
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percentages within a range of 99.4 to 103.8% of the maximum 
net NADH flux were viable. This indicates there are unlimited 
metabolic states that can achieve the FBA identified maximum 
net NADH shunt. The overflow (the percentage above 100%) 
can be offset by the consuming NADH fluxes, for instance, 
the glutamate dehydrogenase (NAD) catalyzed reaction can be 
reversible and consume up to 3.948 (3.757%) of the maximum 
net NADH flux. 

Another significant portion (22.47%) of the maximum net 
NADH flux was supplied by four reactions, catalyzed by alcohol 
dehydrogenase (ethanol to acetaldehyde) (EC: 1.1.1.1), formate 
dehydrogenase (EC: 1.2.1.2), glyceraldehyde-3 -phosphate dehy- 
drogenase (EC: 1.2.1.12), and phosphoglycerate dehydrogenase 
(EC: 1.1.1.95), respectively. These four reactions had rigid flux 
values, which accounted for 5.76%, 7.20%, 5.89%, and 3.62% 
of the maximum net NADH flux respectively. This indicates that 
they are essential to growth. 

All other NADH producing reactions jointly contributed up 
to only 0.06% of the maximum net NADH flux. This suggests 
that, in order to achieve a surplus NADH flux of 105.1 mmol/ 
gDW/h, nearly all of the metabolic resources need to be reallo- 
cated toward the four aforementioned reactions. In addition, the 
rigid variability (v. . /v. > 0.99) of these reaction fluxes under 

<J J * l min i max ' 

high current output indicates that these reactions were much pre- 
ferred for the biomass growth, rather than NADH-targeted cur- 
rent production. 

Furthermore, it was found that the flux through the reaction 
catalyzed by alcohol dehydrogenase (EC: 1.1.1.1) merely had a 
flux of 6.049 mmol/gDW/h under the high degree of NADH 
perturbation. This is contrary to a previous reporting where alco- 
hol dehydrogenase was suggested as the main NADH sink for 
mediators. 15 

Conventionally, glycolysis, TCA cycle and fatty acid oxida- 
tion are considered to be the three main pathways supplying 
NADH to oxidative phosphorylation for ATP generation. 23 
Consequently, it would be appropriate to speculate that these 
three pathways are the potential NADH suppliers for mediators 
in the MFC. However, the presently computed metabolic flux 
model elucidate that the three pathways for NADH regeneration 
are not the metabolic pathways responsible for the excess NADH 
regeneration, subject to the stoichiometry and substrate uptake 
constraints. Instead, two reactions (EC: 1.4.1.3) in nitrogen 
metabolism and (EC: 1.5.1.5) in folate biosynthesis are the only 
two routes allowing a high NADH regeneration rate. 

Aerobic growth. Under aerobic respiration, five reactions were 
identified to have a capability to replenish the loss of cytosolic 
NADH under high current output. Each of the two reactions, 
catalyzed by glutamate dehydrogenase (NAD) (EC: 1.4.1.3) 
and methylene tetrahydrofolate dehydrogenase (NAD) (EC: 
1.5.1.5), had a potential to solely supply up to 81.71% of the 
maximum net NADH flux (113.31 mmol/gDW/h); and any 
combinations of the reaction flux amount to 81.71—81.75% of 
the maximum net NADH flux were attainable. 

About 18.06% of the maximum net NADH flux was supplied 
by alcohol dehydrogenase (EC: 1.1.1.90), glyceraldehyde-3 -phos- 
phate dehydrogenase (EC: 1.2.1.12), and malate dehydrogenase 
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(EC: 1.1.1.85). These three reactions had fixed fluxes, implying 
their enzymatic activities were critical to cell growth. 

Comparison of the identified reactions chosen by anaerobic 
and aerobic growths shows that three reactions, glutamate dehy- 
drogenase (NAD) (EC: 1.4.1.3), glyceraldehyde-3 -phosphate 
dehydrogenase (EC: 1.2.1.12), and methylenetetrahydrofolate 
dehydrogenase (NAD) (EC: 1.5.1.5), were used by the two 
modes to achieve a high NADH regeneration rate. This indicates 
that 87.70% of the maximum net NADH flux in the anaerobic 
mode and 85.50% of the maximum flux in the aerobic mode 
were generated through the same pathways. 

To demonstrate how nutrient uptake is channeled to biomass 
growth and current yield respectively, Table 2 compares the cor- 
responding fluxes at different glucose uptake rates. 

Effect of varying glucose uptake rate on predicted biomass and 
NADH production rates. Table 2 summarizes the linear functions 
obtained for biomass and NADH_mfc production rates vs. glu- 
cose uptake rates. Column 3 shows that the optimal glucose-to- 
biomass efficiency was heavily suppressed under the high current 
output condition in both anaerobic and aerobic growth. 

The linear regression equations are first in units (gDW/ 
gDW/h) (column 4) for comparison with the biomass rate and 
then in the standard flux units (mmol/gDW/h) in column 5. 



Figure 2. The effect of varying biomass production on the fractional 
benefit. The fractional benefit B plotted on the vertical axis, is a mea- 
sure of success in achieving the combined goals of maximal growth rate 
and MET flux. Maximizing one of these at a time, as at the endpoints, 
gives only B = 50%. The graphs show that relative to this, gains in MET 
flux can more than offset losses of growth rate in S. cerevisiae metabo- 
lism. 
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Table 1 . Identified reactions that contribute significantly to the predicted maximum NADH production rate. The flux ranges were calculated from the 
F ATM IN extension of FVA 



Reaction 
ID 


NADH flux 
(mmol/gDW/h) 


Enzyme 


EC No. 


Reaction 


Subsystems 




Min 


Max 










(A) Anaerobic growth 












r_0163 


6.049 


6.049 


alcohol dehydrogenase 
(ethanol to acetalde- 
hyde) 


1.1.1.1 


NAD + ethanol < = > NADH + acetaldehyde 


Pyruvate metabolism 


r_0445 


7.566 


7.566 


formate dehydrogenase 


1.2.1.2 


NAD + formate = > C02 + NADH 


Methane metabolism 


r_0470 


-3.948 


81 .42 


glutamate dehydroge- 
nase (NAD) 


1.4.1.3 


L-glutamate + NAD < = > 2-oxoglutarate + 
NADH + NH3 


Glutamate metabolism 
(aminosugars metabolism) 


r_0486 


6.186 


6.186 


glyceraldehyde-3-phos- 
phate dehydrogenase 


1.2.1.12 


D-glyceraldehyde 3-phosphate + NAD + phos- 
phate < = > 3-phospho-D-glyceroyl phos- 
phate + NADH 


Glycolysis/Gluconeogenesis 


r_0731 


0 


85.37 


methylenetetrahydro- 
folate dehydrogenase 
(NAD) 


1.5.1.5 


5,10-methylenetetrahydrofolate + NAD < = > 
5,10-methenyltetrahydrofolate + NADH 


One carbon pool by folate 


r_0891 


3.808 


3.808 


phosphoglycerate dehy- 
drogenase 


1.1.1.95 


3-phospho-D-glycerate + NAD = > 3-phospho- 
nooxypyruvate + NADH 


Glycine, serine and threonine 
metabolism 


(B) Aerobic growth 


r_0163 


9.971 


9.971 


alcohol dehydrogenase 
(ethanol to acetalde- 
hyde) 


1.1.1.1 


NAD + ethanol < = > NADH + acetaldehyde 


Pyruvate metabolism 


r_0470 


-0.2163 


92.58 


glutamate dehydroge- 
nase (NAD) 


1.4.1.3 


L-glutamate + NAD < = > 2-oxoglutarate + 
NADH + NH3 


Glutamate metabolism 
(aminosugars metabolism) 


r_0486 


4.303 


4.303 


glyceraldehyde-3-phos- 
phate dehydrogenase 


1.2.1.12 


D-glyceraldehyde 3-phosphate + NAD + phos- 
phate < = > 3-phospho-D-glyceroyl phos- 
phate + NADH 


Glycolysis/Gluconeogenesis 


r_0714 


6.190 


6.190 


malate dehydrogenase, 
cytoplasmic 


1.1.1.85 


2-isopropylmalate + NAD < = > 4-methyl- 
2-oxopentanoate + C02 + NADH 


Branched chain amino acid 
metabolism (valine, leucine 
and isoleucine) 


r_0731 


0 


92.58 


methylenetetrahydro- 
folate dehydrogenase 
(NAD) 


1.5.1.5 


5,10-methylenetetrahydrofolate + NAD < = > 
5,10-methenyltetrahydrofolate + NADH 


One carbon pool by folate 



Note: see the supplementary file for the FATMIN results of all the reactions involving NADH and list of the metabolite and reaction abbreviations. 



Comparison of columns 3 and 4 indicates that NADH regenera- 
tion is much less substrate-costly than the biomass production. 

Compared with anaerobic growth, the efficiency of conversion 
of glucose to biomass was 66.67% higher, but the efficiency for 
glucose to NADH was just 7.831% higher in the aerobic growth. 
This implies that the intake of oxygen can promote the biomass 
production but has fewer effects on NADH regeneration. 

Finally, the two operation modes are compared for their theo- 
retically maximum current output in Figure 3 and Table 3. 

Comparison of amperage outputs. Overall, the maximum cur- 
rent output achieved for aerobic growth (6.193 A/gDW) was 
about 7.14% higher than for its anaerobic counterpart (5.781 A/ 
gDW). Since aerobic metabolism can more efficiently use sub- 
strate than the anaerobic, aerobic MFCs were expected to be 
capable of generating more current that the anaerobic one. This 
anticipation is in keeping with the present study's prediction, 
despite the common observation that fuel cells generate less cur- 
rent when they were supplied with substrate (or anolyte) contain- 
ing oxygen. 24 The discrepancy between expectations and reality 



may be ascribed to the deficiency of the engineering design that 
allows the reduced mediators to be intercepted by the presence of 
oxygen rather than anode, resulting in low observed coulombic 
efficiency. Although fermentation produces a lower current than 
oxidative phosphorylation, anaerobic growth eliminates the need 
to aerate large volumes of media and thus is more suitable for 
application in large-scale MFC. 

In general, yeast based MFCs perform better than cyanobac- 
teria but still have a lower power output than bacterial fuel cells. 25 
An experimentally observed current of 16.09 ± 8 uAwas recorded 
when the cell was operating under its lowest impedance in a pre- 
vious S. cerevisiae based MFC. 15 The anode solution contained 2 
g of the yeast, and thus the aforementioned amperage output of 
the MFC system can be converted into 8.45 ± 4 (JtA/g DW, which 
is much lower than the maximum current output (5.714 A/g) 
computed for anaerobic growth in the present study. Another 
experimental study of a methylene blue mediated S. cerevisiae 
MFC reported a maximum current output of 3 mA. 26 With the 
cell density of 50 mg/mL and anode volume of 10 mL provided 
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in the same study, a 0.006A/g can be calculated, 26 which is still 
much lower than the presently computed value (5.714 A/g). The 
large discrepancy between the previously reported current per 
cell gram and the presently computed one may be attributed to 
the fact that (1) the cells were not in optimal metabolic states, 
(2) not every single cell participated in current output, (3) the 
mediator, i.e., methylene blue (MB), could not penetrate the cell 
wall and/or cytoplasmic membrane, and thus could not interact 
with the cytoplasmic metabolic activity efficiently, (4) and the 
engineering design might not efficiently collect all the electrons 
funneled away from the microorganism. 

The eukaryotic respiratory chain has a complicated architec- 
ture which is advantageous to the cell energy balance. 27 The cell 
membrane contains higher protein content, which corresponds 
to the higher rate of redox reactions between cytoplasm and 
inter-membranous space. 27 This feature has been proposed to 
facilitate S. cerevisiae attached to the anode to convey electrons 
to the electric circuit without aids of artificial mediators in an 
MFC. 16 The intracellular redox mediators (i.e., NADVNADH 
and FADVFADH) were suggested to sustain a current output 
of up to 282.83 mA/m. 216 However, without knowledge of the 
cell density on the anode of the MFC, it is impossible convert 
the unit (mA/m 2 ) of the reported current density into a unit of 
amperage per gram dry weight cell, and therefore we are unable 
to compare the current output of the yeast cell in this particular 
MFC to the theoretical maximum value computed here. 

Another more recent study investigated the possibility of 
S. cerevisiae to transfer electrons to an extracellular electron 
acceptor through DET mode and found that the cells adhered 
to the anode were able to sustain power generation in a media- 
tor-less MFC configuration. It is proposed that the surface con- 
fined species are responsible for the DET, 14 however, the power 
performance of this MFC was extremely low (0.003 W m~ 2 and 
0.03 A/m 2 ). 14 With a reported cell sample weight of 0.13 g and 
the anode surface area of 4 cm 2 , the current density can be trans- 
lated into 0.000009231A/g. Since the DET mode is based on 
surface confined species and it is unclear what those species are, 
identification of the protein involved in such a DET mode is 
required, before the in silico method presented in this study can 
be employed to reveal the innate capability of the DET mode 
based on S. cerevisiae. 

Conclusions 

In this study, the anaerobic and aerobic growth modes of S. cere- 
visiae were modeled to study its maximum potential for MFC 
current output and underlying metabolic metabolisms. It was 
shown that a similar level of maximum current outputs can be 
achieved for the two growth modes (i.e., 5.781 A/gDW for the 
anaerobic growth and 6.193 A/gDW for the aerobic growth), 
whereas under aerobic conditions the yeast achieved nearly twice 
the anaerobic growth rate. To achieve the highest efficiency of 
NADH regeneration and maximize the current output, the aero- 
bic and anaerobic metabolisms relied on two reactions, catalyzed 
by glutamate dehydrogenase (NAD) (EC 1.4.1.3) and methylene 
tetrahydrofolate dehydrogenase (NAD) (EC 1.5.1.5), for NADH 



7 




Electron flux (mmol gDW 1 h 1 ) 



production. Besides, the robustness analysis indicated that the 
NADH regeneration was much less energy costly than the bio- 
mass production rate. Taken together, our findings suggest that 
S. cerevisiae should be re-evaluated and receive more research 
effort for MFC electricity generation. 

Materials and Methods 

General modeling assumptions. The modeling was conducted 
based on the three assumed conditions: (1) The MFC reac- 
tor is a chemostat that can provide all the necessary nutrients 
to maintain the optimal microbial growth; (2) Only a pure 
S. cerevisiae culture is used for electricity generation; (3) A puta- 
tive ideal mediator (e.g., Neutral red 28 ) can be found to enter the 
cytoplasm, intercept the electrons from NADH-involved meta- 
bolic activity and convey the electrons to the anode; and (4) The 
uptake rate of the organic substrate (i.e., glucose) is allowed to 
vary but is constrained to a realistic range to reflect the observed 
uptake rates for optimal growth. 

Modeling electrode interactions. The interactions with an 
electrode were captured by introducing two reactions into the 
model reconstruction (Table 4) . These reactions represent the net 
reaction between the reducing equivalents and the electrodes in 
the MFC. 

The two added reactions reflect a well-established understand- 
ing 29 that mediators in an MFC enter the cell and interact with 
cytosolic NADH in the metabolism to extract energy (in the 
form of electrons) from these electron shuttles, as schematically 
shown in Figure 4. Because NADH is a major electron carrier 
in microbial metabolism connecting many different cellular pro- 
cesses, known as a "hub" in terms of network theory, 30 the effects 
of the NADH drain that the added reactions creates, propagates 
throughout the metabolic network. In this way, the modeling can 
estimate the maximum power achievable for a microorganism in 



Figure 3. The current output (A/g) as a function of electron flux. The 
dark red line denotes the maximal current outputs and NADH_mfc 
production rate, while the area represents all allowable current outputs 
and electron production rates. The round dotted arrow line indicates 
the maximal current output and corresponding electron production 
rate when the growth rate is set to 5% of the predicted maximum 
growth rate (0.01303 h" 1 ). 
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Table 2. The summary of linear functions of the biomass production and MET rates in the heterotrophic, photoautotrophic, and mixotrophic modes 



Metabolic type 


Substrate 


Growth 




MET 






(gDW/gDW/h) 


(gDW/gDW/h) 


(mmol/gDW/h) 


Anaerobic 


glucose 


y c = 0.0261 x 


y, , „„„,, = 6.971x 

' (g/gDW/h) 


y, „ n,.„ k , = 10.509X 

' (mmol/gDW/h) 






y =0.0035x 
'p 






Aerobic 




y c = 0.0435 x 


y, , „„„.. = 7.518x 

7 (g/gDW/h) 


y. „ = 11.332x 

' (mmol/gDW/h) 






y =0.0038x 







Note: The FBA simulations were performed by changing the photon uptake rate and the glucose uptake rate with the maximization of the objectives, 
y denotes the MET rate (mmol/gDW/h or gDW/gDW/h) or growth rate (gDW/gDW/h), whereas x represents the substrate (i.e., glucose) uptake rates 
(mmol/gDW/h). y c represents the optimal growth, whereas y p denotes The perturbed growth under cytosolic NADH deprivation. 



Table 3. Comparison of predicted amperage output under theoretical maximum current output condition and 5% of optimal growth rate condition 

Mode Conditions Biomass production rate Electron (mmol gDW 1 h ') Amperage (A gDW 1 ) Coulombic efficiency (CE%) WgDW 1 

(h 1 ) 

Anaerobic 5% of optimal 0.01303 215.7 5.781 89.87% 4.798 

Aerobic growth rate 231/| g _ lg3 96.29% 5.141 



Electric circuit 
< — 




. Cell wall 
Plasma membrane 



Mitochondrion 



Figure 4. A schematic of the modeled anodic mechanisms where the cytoplasmic NAD7NADH 
cycle as the electron supplier. Microbes take up substrates (glucose) generating carbon dioxide 
and proton. This process yield electrons for metabolic benefit, i.e., growth, and reduces Med ox in 
the cytosol into Med ,. Med , diffuses into contact with the electrode, where Med .reduces the 

' red red ' red 

electrode generating electrical current. The oxidized form, Med , diffuses back through anolyte 
for reuse by the microbes. The radial red circle highlights the redox cycle supplying the electron to 
the anode, which is investigated in the present modeling. 



maximization of the production of 
an overflow flux of desirable reducing 
equivalent, cytoplasmic NADH, so as 
to restore the redox balance perturbed 
by the energy extraction process during 
MFC operation. The bi-level objective 
optimization can result in a set of Pareto- 
optimal points, 37 forming a Pareto front, 
which describes the trade-off relationship 
between the growth rate of the organism 
and the regeneration rate of a plethora 
NADH flux. Through varying the rela- 
tive weights associated with each objec- 
tive, 38 the impact of the enhanced energy 
extraction on the growth of the cell can 
be evaluated. The bi-objective equation 
in this study is defined as follows in terms 
of the biomass growth flux F B and the 
reducing equivalent diversion flux F N : 

0=(l-\)F B + \gF N 

A unit conversion factor g is explic- 
itly provided because in the model, F N 



MFC, as well as elucidate the response of energy metabolism to 
the perturbation induced by electricity generation. 

Objective equation. The objective equation was formulated 
as described in our recent studies. 19,22 Briefly, we constructed a 
bi-objective equation that comprises two objectives: one is the 
maximization of the growth rate, which reflects the belief that 
microbes evolve to optimize their metabolisms for maximiz- 
ing biomass production (growth rate) 31 " 34 and has been experi- 
mentally validated by numerous studies 32,35,36 ; the other is the 



and F B are measured respectively in units 
(mmol/gDW/h and g/gDW/h or [h 1 ]). 
This ensures that X is a dimensionless 
fraction that can be directly interpreted 
as the relative contribution of the NADH 
flux to the combined objective. 

Fractional benefit analysis. The bi-objective optimization 
conducted in this study can assess how metabolic energy loss dur- 
ing current production influences the biomass growth. However, 
the mathematical optimization method based on such a bi-objec- 
tive function cannot clarify the combined benefit resulting from 
an optimized metabolic state. Therefore, this requires implemen- 
tation of fractional benefit analysis to mathematically quantify 
either extreme on the Pareto front; even though only one of the 
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objective terms is maximized, the other term also has a value and 
will contribute (positively or negatively) to the overall benefit. 
In the fractional benefit analysis, each of F B and F N is scaled as a 
fraction of the maximal values F B * and F N * they can ever achieve, 
i.e., the values when each is independently maximized. As the 
maximum of each of these fractions is 1, a plausible measure of 
the combined benefit achieved in any particular metabolic state 
is the value of the quantity we denote as the fractional benefit B: 



B = 



i 7 ,, F„ 



-+- 



A value of B in excess of 50%, indicates that gains in one 
objective more than compensates for losses in the competing 
objective. 

Conversion of units of flux and current. Current (in amperes) 
was integrated over time and converted to electrons recovered by 
using the following conversions: 1 C = 1 A x 1 s, 1C = 6.24 x 10 18 
electrons, and 1 mol = 6.02 x 10 23 electrons (the Faraday constant 
96485 C/mol). Therefore, one flux unit (mmol/g/h) can be con- 
verted into A/g as follows: 

1 1/ /u lmo1 96485 C nn^o a / 

1 mmol/g/h = x = 0.0268 A/g 

lOOOgx 3600s mol 

Calculation of Coulombic efficiency (CE). One of the 

parameters commonly used to quantify the performance of 
MFCs is the Coulombic efficiency (CE). The CE is defined as 
the ratio of electrons transferred to the anode to that in the start- 
ing substrate. We use the full oxidation of glucose with oxygen 
as the oxidant as the reference reaction to characterize the energy 
efficiency of the respiratory metabolism: 

Glucose oxidation reaction: C,H..O, + 6 FLO -> 6 CO^+ 24 
FT+24 e" 

Based on such stoichiometric information, 1 mol of glucose 
supplies 24 mol of electrons. 39 



CE%= °"' p '" xl00% = 
C 



The MET flux(mmol/gDW/h)x 100 c 
acetate uptake rate (mmol/gDW/h) x 24 



Calculation of theoretical power outputs of the three tested 
modes. An upper limit for the cell voltage is calculated in this 
work based on formal potentials of the biological and electro- 
chemical redox processes, as given by: 



AE 



cell 



cathode 



anode 



AE 



Here, cel1 is the standard cell potential (aka., electromo- 

E°' 

tive force); cathode j s tne standard potential of cathode oxida- 

E°' 

tion; anode is the standard potential of anode reduction. The 
formal potentials of the anode and cathode used for calculation 
of power density in the three operation modes are summarized 
in Table 5. 

The real potential derived from the MFC will be reduced 
due to various potential losses 12 associated with MFC operation, 
such as ohmic resistances, concentration polarization and kinetic 



Table 4. The added reactions for modeling the interaction of the micro- 
organism and the electrode in MFCs 



Operation mode of the Reaction ID 


Reaction 


MFC 





MET INADHmfc nadh -> nadh_mfc 

2NADHmfc nadh_mfc -> nad + h_emm 

nadh_mfc: the NADH available for MET mode of MFCs; 

h_emm: the H ions as the by-product released from the reaction 
2NADHmfc; 



Table 5. The standard potential of the redox reactions involving the 
electron donor and acceptor for the MFC (measured at pH 7). 







Redox couple 


E°'(V) 


Anode 


MET 


NAD + +H + +2e^NADH 


-0.320 47 


Cathode 




0 2 +4H + +4e^2H 2 0 


+a51 48-50 



Table 6. The theoretical limit of standard anode potentials of MFC based 
on S. cerevisiae 



Electron transfer 
mode 


"^anode 


^cathode 




MET 


-0.32 


0.51 


0.83 



constraints. These are not taken into account in the calculations 
reported here. 

The MFC standard cell potential calculated as above and 
shown in (Table 6) shows that as long as the same electron donors 
are used, the choice of biocatalysts will have little effect in the cell 
potential. 

Calculation of the viable growth rate. To obtain a realistic 
upper value for the current that can be extracted from a microbe, 
we chose 5% of the maximum theoretical biomass production 
rate as the minimum viable growth rate. This percentage value is 
more likely to be achievable in practice, when compared with 1% 
assumed as the viability threshold for computational identifica- 
tion of the lowest growth rate in several previous studies. 40,41 

Simulating S. serevisiae growth. A recently published met- 
abolic network of S. cerevisieae, Yeast 5 GEM 21 (Version 5.40; 
November 2012), was chosen as the analysis backbone for the 
all the simulations. Compared with previous metabolic networks 
of this yeast, this network model is a fully compartmentalized 
and elementally-balanced S. cerevisiae metabolic network, which 
includes more genes and reactions based on genomic, biochemical 
and physiological information. The details of existing reactions 
(substrate and cofactor specificity, reaction reversibility, and com- 
partmentalization) in the network have also been re-evaluated to 
update the model based on existing literature. Furthermore, the 
SBML file compatible for computation with the COBRA toolbox 
is attached to the publication. 42 

All simulations were performed on glucose as the sole carbon 
source. The exchange of ammonium, protons, iron(2*), phos- 
phate, potassium, sodium, sulfate, and water were all uncon- 
strained. The maximum rate of the glucose uptake is set to 



428 



Bioengineered 



Volume 4 Issue 6 



10 mmol/gDW/h. 43 For aerobic growth, oxygen exchange was 
unconstrained. 21 For simulation of anaerobic growth, the oxygen 
exchange reaction was constrained to 0. In addition, ergosterol, 
lanosterol, zymosterol, and phosphatidate were allowed for free 
exchange and the biomass definition was modified by removing 
14-demethyllanosterol and ergosta-5,7,22,24(28)-tetraen-3beta- 
ol from the "lipid" definition. 21 These steps were applied to reflect 
the observation that yeasts require sterols and fatty acids when 
cultured under rigidly anaerobic conditions. 21 

Analysis technique. The growth rates and electron pro- 
duction were computationally determined using Flux Balance 
Analysis (FBA). 44,45 Computations were performed with 
COBRA Toolbox 42 in MATLAB (The Math-Works Inc) and 
OptFlux. 46 Flux variability analysis with target flux minimiza- 
tion (FATMIN) 20 was conducted to determine the flux ranges 
of the reactions in the network under the metabolic states heav- 
ily perturbed by electricity generation. FATMIN can examine 
the feasible flux ranges of the desired reactions (such as NADH- 
involved ones) of a genome-scale metabolic network, while elimi- 
nating the pointlessly high fluxes values of the loop reactions, to 
achieve biologically meaningful values. 

The underlying idea of FATMIN is to collect the fluxes of all 
reactions that produce the target metabolite (as opposed to those 



that consume it) and minimize this collective flux while leaving 
the value of the objective function unchanged. In this way, the 
flux in reactions that contribute both productive and futile flux 
components is limited to just the productive part. Because of the 
flux balance conditions, futile fluxes are eliminated also from all 
other reactions that form part of the cycle. 

FATMIN is a pipelined algorithm comprising FBA and flux 
variability analysis (FVA) 20 and therefore can be performed with 
either COBRA or OptFlux according to users' preference. The 
maximum reducing equivalent production rates for arbitrary 
growth rates between the data points were calculated using the 
interpolation function in Mathematica 8.0 (Wolfram Research, 
Inc). All the FBA and FATMIN results are detailed in the 
Supplemental Tables. 
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